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摘 要 : DOME SSN EWA UR MAAR AY FIERA, SEBS RMU EKRA a AURA 
VOTE AA BH BU Ae EI He IZA AA A PA, E 27.6*C， 相 对 湿度 34.7% 空 气 入 口 条 件 下 分 别 模拟 了 在 ISP. Bh 
液 层 及 组 织 导 热 、 对 法 和 汽化 潜 热 和 2) 呼吸 道 壁面 恒定 37°C 两 种 边界 条 件 下 的 理想 et 颗粒 物 的 运动 和 
沉积 。 结 果 表 明 ， 条 件 1 平均 相对 湿度 较 条 件 2 高 2.3%， 但 笑 粒 物 沉积 率 比 条 件 2 小 ， 最 大 差异 可 达 9%; 颗粒 物 沉 积 形 式 类 
似 ， 但 条 件 1 下 沉积 的 颗粒 物 平 均 直 径 比 条 件 2 小 40.7%， 逃 逸 颗粒 物 平均 直径 小 14.3%。 
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Numerical study of the transport and deposition of hygroscopic particles 


in a mouth-throat airway with heat-transfer boundary 
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Abstract: 

A mouth-throat (MT) airway model is constructed, including the airway tissue and mucus layer. It realizes the 
simulations of the water vapor transport above the mucus layer, calculations of the convection heat-transfer, the 
latent heat due to evaporation of the mucus and the heat condition within the airway tissue. The transport and 
deposition of multi-component hygroscopic particles are simulated at 27.6°C and relative humidity (RH) at 34.7%. 
Two types boundary conditions are considered, 1.e., 1) considering heat transfer, including the conduction, 
convection and latent heat generation in the airway, mucus layer and tissue, 2) constant temperature of the airway 
boundary at 37 °C. The simulation results indicate that the average RH of condition 1 is 2.3% higher when compared 
to that of condition 2. However, the deposition efficiency (DE) of condition 1 is smaller, and the largest DE 
difference between the two conditions is up to 9%. The deposition patterns are similar when comparing the two 
conditions. But the average diameters of deposited and escaped particles of condition 1 are 40.7% and 14.3% smaller 
than that of condition 2, respectively. 
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058 长 ， 直 径 和 惯性 发 生 明显 提高 ， 使 颗粒 物 能 够 高 效 


近年 来 随 着 呼吸 道内 可 吸入 颗粒 物 研究 的 不 
断 深入 ， 研 究 者 越 来 越 关 注 利 用 可 吸入 颗粒 物 运动 
和 沉积 规律 尝试 实现 呼吸 道内 药物 颗粒 的 肢 问 投 
递 ， 其 中 吸湿 性 颗粒 物 相 关 研 究 引 起 了 较为 广泛 的 
关注 。 已 有 研究 表明 ， 通 过 控制 吸湿 性 颗粒 物 的 初 
始 直 径 ， 使 其 在 罕 过 上 呼吸 道 时 较 少 地 发 生 沉积 ， 
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地 在 肺 内 沉积 。 

由 于 吸湿 性 甘 粒 物 的 吸湿 增长 过 程 多 受 坏 境 
温度 和 湿度 影响 了， 所 以 准确 模拟 呼吸 道内 温度 、 
湿度 条 件 是 预测 吸湿 性 可 吸入 箱 粒 物 运动 和 沉积 的 

前 提 条 件 。 然 而 ， 在 以 往 的 研究 中 一 般 使 用 呼吸 道 
璧 面 温度 为 37?C、 相 对 湿度 (Relative Humidity, RH) 
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为 99.5% 或 100% 条 件 进行 研究 下 。 同 时 目前 的 实 
验 中 一 般 采 用 3D 打印 和 快速 成 型 方法 加 工 呼 吸 道 
模型 ， 材 料 一 般 为 塑料 、 树 脂 等 ， 因 此 难以 模仿 呼 
吸 道 表面 香 液 输 运 和 组 织 传 热 过 程 * 小, 但 实际 过 程 
中 ， 周 围 环 境 温 度 和 RH 一 般 低 于 呼吸 道 条 件 ， 所 
以 在 呼吸 过 程 中 ， 呼 吸 道 壁 面 黏 液 层 及 附近 组 织 将 
与 空气 进行 传 热 ， 包 括 与 空气 的 对 流 换 热 、 水 分 蒸 
发 所 造成 的 气 化 潜 热 损失 和 组 织 内 导热 ”。 

本 文 构建 了 含有 科 液 屋 和 周边 人 体 组 织 的 理 
想 口 喉 模 型 ， 并 在 不 同 传 热 边界 条 件 下 模拟 了 吸湿 
性 可 吸入 颗粒 物 的 运动 和 沉积 过 程 ， 即 在 空气 7 = 
27.6"C，R 琳 = 34.7% 入 口 条 件 下 对 比 了 1) 考虑 组 织 
内 导热 、 符 液 层 表面 对 流 和 水 蒸气 蒸发 造成 的 汽化 
潜 热 损 失 以 及 2) 呼 吸 道 壁 面 恒定 37°C 两 种 边界 条 
件 下 的 温度 和 相对 湿度 分 布 以 及 多 组 分 颗粒 物 的 沉 
积 率 和 沉积 形式 。 
1 控制 方程 
1.1 气相 

由 于 口 喉 模型 中 存在 横 规 面 渐 缩 ， 流 速 逐渐 增 
大 ,存在 层 流 - 消 流 换 变 和 消 流 流动 ,因此 气相 场 流 
动 模拟 采用 了 Transition SST 模型 ， 因 篇 幅 限 制 不 
罗列 控制 方程 ， 细 节 可 参考 前 人 研究 下 ， 在 呼吸 
道内 的 应 用 可 参考 Zhang 和 Kleinstreuer(2011) 工 作 
中。 为 分 析 水 共 气 在 呼吸 道内 的 分 布 ， 采 用 了 干 空 
气 和 水 蒸气 两 种 组 分 ， 组 分 输 运 模型 为 -: 


<(pY.)+V-(p0¥,) =-V-J ( 


相应 能 量 方程 为 “: 


A ER =V [Eva] (2) 
Ot C 
1. 2 颗粒 相 

颗粒 相 控制 方程 主要 包含 两 方面 ， 一 方面 是 运 
STi Re, ABCC ATR ACUI Pa, RR RUAN TY tal 
撞 ， 其 运动 可 由 牛顿 第 二 定律 确定 “|; 


mdi, /dt = F, +m,(p, — p)z/ P, (3) 


Fy Fy Ee A E a LTS 28 A EE 
方程 ， 该 模型 假设 颗粒 物 表 面 有 一 厚度 可 忽略 的 注 
屋 ， 湾 层 中 气态 组 分 浓度 由 轿 粒 物 ( 液 滴 ) 表 面 组 分 
决定 。 湾 层 内 又 气 成 分 与 周围 流体 中 北 气 成 分 的 传 


质 过 程 中 ,可 燕 发 组 分 e 的 平均 质量 通 量 为 U191, 


(4) 


由 此 可 以 确定 吸湿 性 颗粒 物 吸 湿 增 长 /蒸发 "; 


dm, k k 二 
P > | 1.d4 x p2 (ne - A) (5) 
AOFERE e WREDA Y, ur PEAH 
修正 的 Raoult 定律 得 到 1> 14, 
Y — K l (T,) 6 
e,surf = Toe ”DRT ( ) 
1.3 固 相 
固 相 中 仅 为 热传导 过 程 ， 其 控制 方程 为 : 
<(phi)=V-(kVT) +5, (7) 
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2 模拟 对 象 与 条 件 


理想 口 喉 模 型 空 及 部 分 结构 与 Zhang 等 
(2006 1 所 提出 的 相同 ， 口 腔 入 口 直径 为 30mm, 
喉 部 出 口 为 8.5mm。Wnu 等 (2014) 研 究 指 出 , 受 传 热 
影响 的 呼吸 道 周围 组 织 厚度 约 0.5mm", AEE 
Zhang 等 (2006) 所 用 模型 外 增加 了 总 厚度 为 1mm 的 
蔚 液 层 -组 织 区 域 ,理想 口 喉 模型 结构 和 网 格 划 分 示 
意图 如 图 1 所 示 ， 其 中 流体 域 采 用 了 Ogrid， 进 行 
了 边界 层 加 密 ， 保 证 首 层 网 格 y 值 小 于 1。 划 分 和 
液 层 -组 织 区 域 网 格 时 ， 采 用 了 类 似 的 加 密 方法 ， 贴 
近 流体 的 第 一 层 网 格 厚 度 为 10um, ERRE, 
其 余 九 层 网 格 逐 渐 增 厚 ， 占 据 990um， 作 为 呼吸 道 
周围 组 织 。 生 成 了 网 格 数 为 61.6 万、289.3 万 和 988.5 
万 的 三 个 网 格 ， 经 网 格 无 关 性 校 验 ， 选 择 网 格 数 为 
289.3 万 的 网 格 ， 其 中 流体 域 和 黏液 层 -组 织 区 域 网 


格 数 分 别 为 228.7 万 和 60.6 万 。 

在 模拟 中 假设 稳 态 吸 气 流量 为 23.6 L/min( 入 口 
处 雷诺 数 Re = 832.6， 出 口 处 Re = 2938.7)， 入 口 处 
速度 分 布 为 层 流 的 抛物 面 形状 ， 假 设 壁 面 RH = 
99.5%， 入 口 空气 T= 27.6"C，R 厅 = 34.7%， 首 先 在 
1) 考 上 处 导热 、 对 流 和 汽化 潜 热 和 2) 呼吸 道 壁 面 恒定 
37°C 两 种 条 件 下 开展 干 空气 和 水 蒸气 两 种 组 分 气 
体 的 流动 模拟 。 


流 回 界面 


图 1 理想 口 喉 模 型 与 网 格 划分 示意 图 
Fig. 1 Idealized MT airway model and mesh 


在 流动 模拟 收敛 后 ( 残 差 < 10") 在 稳 态 流 场 下 
进行 吸湿 性 多 组 分 颗粒 物 运动 和 沉积 的 模拟 ， 颗 粒 
物 初 始 组 分 为 水 、 酒 精 、NaCl 和 荧光 剂 ， 质 量 比 为 
400: 100: 100: 2.5， 在 低 RH 条 件 下 ， 颗 粒 物 将 
蒸发 失去 水 和 酒精 成 分 ， 形 成 固体 颗粒 ( 仅 含 NaCl 
和 荧光 剂 ， 直 径 为 初始 直径 的 44.3%)， 当 颗粒 物 进 
入 壁面 附近 的 高 R 区 域 时 可 能 发 生 吸湿 增长 。 如 
表 1 所 示 ， 每 种 边界 条 件 进 行 8 组 不 同 颗粒 物 初 始 
直径 的 模拟 以 获得 颗粒 物 沉 积 率 与 Stokes 数 关 系 。 
每 次 颗粒 物 运 动 和 沉积 模拟 中 在 入 口 处 释放 10000 
个 颗粒 ， 颗 粒 物 在 入 口 的 分 布 的 概率 密度 正比 于 当 
地 速度 ， 即 为 抛物 面 型 ， 颗 粒 物 初 始 速度 与 流速 相 
同 。 为 保证 颗粒 物 吸湿 /蒸发 模拟 的 稳定 ， 颗 粒 物 运 
动 的 时 间 步 长 设 定 为 10s。 此 外 ， 颗 粒 物 的 脉动 速 
度 “和 灌流 脉动 的 近 壁 面 修正 "由 自 编 UDF 完成 。 


R1 吸湿 性 多 组 分 颗粒 物 初始 直径 与 Stokes 工 况 表 
Table 1 Initial diameters and Stokes numbers of the 


multi-component hygroscopic particle 


颗粒物 直径 d, (um) Stokes 数 
二 1.98E-04 
3 1.24E-03 
7.5 2.79E-03 
10 4.96E-03 
11.25 1.12E-02 
12.5 1.98E-02 
13.75 2.51E-02 
15 3.10E-02 


3 RAW 


3.1 温度 与 RH 分 布 
由 于 篇 幅 限 制 ， 在 此 不 详细 介绍 理想 口 喉 模型 

内 的 流 场 分 布 ， 细 节 可 参见 前 人 研究 11。 在 图 2 给 
出 传 热 边界 下 (条 件 1) 和 恒定 壁 温 ( 条 件 2) 中 平面 、 

出 口 以 及 黏液 层 中 的 慢 度 分 布 。 从 整体 中 平面 和 出 
Sa Cte wears 
分 相似 ， 但 值 的 大 小 相反 ， 这 是 由 于 进口 空气 温 
低 于 呼吸 道 边界 ， 千 We 
地 与 呼吸 道 表 面 发 生 热 交换 。 考 虑 呼吸 道 表 面 导 热 、 
对 流 和 汽化 潜 热 作用 时 (图 2(b))， 能 够 明显 观察 到 
秋 液 层 中 温度 的 不 均匀 分 布 ， 在 口 咽 附 近 的 这 区 处 
ee 蒸发 ， 因 此 
温度 较 低 。 在 恒定 壁 温 条 件 下 ， 可 以 观 聚 到 壁面 附 
oe H E OAE RAE- 
组 织 内 并 无 温度 梯度 。 
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(b) 条 件 1 黏液 层 温 度 分 布 
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(c) 条 件 2 中 平面 及 出 口 温度 分 布 
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(d) 条 件 2 ABE in Ea Ah 


2 边界 传 热 和 恒定 壁 温 条 件 下 口 喉 模型 内 中 的 温度 分 布 
Fig. 2 The temperature distributions of the MT airway for 
heat-transfer boundary and constant temperature boundary 

conditions 
图 3 所 示 为 两 种 条 件 下 理想 口 喉 模型 内 中 平面 
及 出 口 的 RH 分 布 ， 两 者 较为 接近 ， 但 能 够 友 现在 
口 咽 处 和 出 口 处 ， 条 件 1 下 的 高 RH 区 域 更 深入 地 
扩展 到 窒 流 中 部 区 域 ， 对 比 两 者 在 流体 域内 的 平均 
RH 可 以 及 现 条 件 1 比 条 件 2 高 2.3%。 在 前 期 研究 


中 已经 友 现 ， 吸 湿性 壬 粒 物 的 增长 过 程 对 RH 较为 
敏感 "， 如 果 扩 展 到 整个 呼吸 系统 ， 可 能 带 来 较 大 
的 吸 温 性 颗粒 物 沉 积 差 寞 。 


= 


(a) 入 口 空气 7= 27.6°C，RH= 34.7%， 考 虑 导热 、 对 流 和 
汽化 潜 热 ， 流 体 域 平均 RH = 56.9% 


(b) 入 口 空气 7=27.6°C, RH = 34.7%, 呼吸 道 壁面 恒定 37°C, 
流体 域 平均 RH = 54.6% 

图 3 不 同 条 件 下 理想 口 喉 模型 内 中 平面 及 出 口 的 RH 分 布 

Fig. 3 RH distributions of the symmetrical plane and outlet of 


the idealized MT model for different boundary conditions 

3.2 沉积 率 

图 4 所 示 为 两 种 条 件 下 吸湿 性 颗粒 物 沉 积 率 与 
Stokes 数 关 系 ， 可 以 上 友 现 在 条 件 1 下 ， 流 体 域 平均 
RH 相对 较 噩 ( 见 图 3(a)， 但 其 沉积 率 仍 低 于 条 件 2 
的 结果 (Sf=0.025 时 ， 沉 积 率 差异 最 大 约 为 9%)。 这 
与 以 往 在 定 温 条 件 下 颗粒 物 吸湿 能 力 与 RH 正 相关 
有 所 差异 , 说 明 吸 湿性 壬 粒 物 吸湿 增长 / 芭 友 过 程 可 
能 不 仅 与 RH 有 关 ， 还 与 温度 有 关 。 因 此 ， 人 简化 恒 
定 壁 温 条件 可 能 对 高 估 吸 湿性 颗粒 物 沉积 京 ， 同 时 
在 此 条 件 下 呼吸 道内 大 部 分 区 域 RH 疝 低 于 NaCl 
吸湿 增长 的 阔 值 …”， 肺 内 高 RH 环境 的 模拟 仍 有 
EF JE o 
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图 4 不 同 条 件 下 吸湿 性 可 吸入 颗粒 物 沉积 紊 与 Stokes BK 
系 


Fig. 4 Comparison of deposition efficiency of hygroscopic 


particles under different boundary conditions 


3. 3 沉积 分 布 

图 5 所 示 为 两 种 条 件 下 初始 直径 11.2Sum 吸湿 
性 颗粒 物 的 沉积 和 离开 出 口 时 的 分 布 及 最 终 颗 粒 直 
径 。 由 于 两 种 条 件 的 差异 主要 在 温度 上 ， 所 以 流 场 
以 及 颗粒 物 轨迹 基本 相近 。 可 以 发 现 ， 沉 积 的 颗粒 
物 大 多 集中 在 咽喉 区 域 的 这 曲 处 ， 同 时 颗粒 物 相对 
均匀 的 沉积 在 了 咽喉 一 周 ， 说 明 颗 粒 物 一 定 程度 上 
跟随 了 二 次 流 的 运动 ， 而 不 是 直接 受 惯性 控制 碰撞 
在 后 侧 壁 面 ， 这 应 与 颗粒 物 在 入 口 段 低 RH 区 域 水 
和 酒精 组 分 蒸发 导致 颗粒 变 小 有 关 。 同 时 观察 出 口 
处 的 颗粒 物 分 布 , 发现 大 笑 粒 大 多 集中 在 壁面 附近 ， 
即 高 RH 区 域 ， 部 分 大 颗粒 被 二 次 流 卷 吸入 气管 中 
部 , 其 余 RE<70%( 人 参考 图 3) 区 域内 颗粒 物 基 本 为 干 
颗粒 (图 5 FERE). Sita 5 中 的 颗粒 物 直径 得 
知 条 件 1 下 沉积 的 颗粒 物 平 均 直 径 为 6.20pm, 逃逸 
的 颗粒 物 平 均 直 径 为 6.15um, fF 2 下 分 别 为 
10.45um 和 7.18um。 这 种 差异 可 能 在 沉 场 复 林 的 蜡 
腔 或 更 大 的 范围 内 的 呼吸 道内 变 得 更 加 明显 。 
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(a) 考虑 导热 、 对 流 和 汽化 潜 热 


e 10 um 


(b) 呼吸 道 壁面 恒定 37°C 
图 5 不 同 条 件 下 理想 口 喉 模 型 内 初始 4,= 11.25um 吸湿 性 
颗粒 物 沉积 形式 和 出 口 分 布 


Fig. 5 Deposition patterns and distributions at the outlet of 


hygroscopic droplets with initial diameter of 11.25um in the 


idealized MT model under different boundary conditions 
十 
4 2 论 


本 文 构 建 了 含有 黏液 屋 和 周边 人 体 组 织 的 理想 
口 喉 模 型 ， 在 27.6*C， 相 对 温度 34.7% 空 气 入 口 条 
件 下 分 别 模拟 了 在 1) 考 虑 呼吸 道 、 医 液 层 及 组 织 导 
热 、 对 流 和 汽化 潜 热 和 2) 呼吸 道 壁 面 恒定 37°C 两 
种 边界 条 件 下 开展 了 多 组 分 吸湿 性 可 吸入 颗粒 物 的 
运动 和 沉积 模拟 ， 对 比分 析 了 理想 口 喉 模型 内 的 温 
EM RH 分 布 以 及 颗粒 物 的 沉积 率 和 沉积 形式 ， 主 
要 结论 如 下 : 

1) 当 考虑 了 导热 、 对 流 和 汽化 潜 热 后 ， 能 够 明 
we LSE BIA A i PEK 37°C 和 温度 的 不 均匀 


分 布 ; 3) 两 种 条 件 下 颗粒 物 沉 积 形 陈 相似 ， 但 条 件 2 
2) 条 件 1 下 平均 RH 较 条 件 2 ZAI 2%, 但 颗粒 下 沉积 和 逃逸 的 颗粒 物 的 平均 直径 比 条 件 1 分 别 高 


物 沉积 率 反 而 较 低 ， 最 大 差异 可 达 9%， 说 明 颗 粒 40.7% 和 14.3%， 这 一 差异 在 复杂 流 场 (如 蜡 腔 ) 和 大 
物 的 吸湿 增长 不 仅 与 RH FAK, 还 和 温度 分 布 有 关 ; YE FPR HEB he NBA, WA FEIT RE 
主要 符号 表 

Nomenclature 

C, 比 定 压 热 容 m, 颗粒 质量 

d, 颗粒 直径 Pa(T,) ”7 下 的 饱和 压力 

Di, 水 蒸气 热 扩 散 系数 Sc, Schmidt 数 

D. 组 分 e 的 质量 扩散 率 Sh Sherwood 数 

E, 颗粒 电力 u, 颗粒 速度 

H SIKIA x, 组 分 e 在 溶液 中 的 摩尔 分 数 

J P ROEE Y 气体 质量 分 数 

k 导热 系数 Ve 组 分 e 在 溶液 中 的 活 度 系数 
k, a VL AG EE 4, WL JIE 

K, 组 分 e 的 Kelvin 效 应 修正 系数 P, POUL PE 
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